Anomalous Dynamics of Translocation 
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We study the dynamics of tfie passage of a polymer through a membrane pore (translocation), 
focusing on the scaling properties with the number of monomers A'^. The natural coordinate for 
translocation is the number of monomers on one side of the hole at a given time. Commonly used 
models which assume Brownian dynamics for this variable predict a mean (unforced) passage time 
T that scales as A''^, even in the presence of an entropic barrier. However, the time it takes for a free 
polymer to diffuse a distance of the order of its radius by Rouse dynamics scales with an exponent 
larger than 2, and this should provide a lower bound to the translocation time. To resolve this 
discrepancy, we perform numerical simulations with Rouse dynamics for both phantom (in space 
dimensions d = 1 and 2), and self-avoiding (in d = 2) chains. The results indicate that for large 
A'^, translocation times scale in the same manner as diffusion times, but with a larger prefactor that 
depends on the size of the hole. Such scaling implies anomalous dynamics for the translocation 
process. In particular, the fluctuations in the monomer number at the hole are predicted to be 
non-diffusive at short times, while the average pulling velocity of the polymer in the presence of a 
chemical potential difference is predicted to depend on A'^. 
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I. INTRODUCTION 

The process of translocation, in which a polymer 
worms its way through a narrow pore, is an event im- 
portant to many biological systems. Examples include 
the viral injection of DNA into a host, DNA packing into 
a shell during viral replication, gene swapping through 
bacterial pili, and the genetic technique of cell transfor- 
mation by DNA electroporation There are also a 
number of recent in vitro experiments on translocation, 
e.g. the electric field-induced migration of DNA through 
microfabricated channels |^ or through an a-hemolysin 
protein channel in a membrane |^,^. The driving force 
is an essential ingredient in the above process, as are 
the entropic and cooperative factors that arise from the 
connectivity of the polymer. An interesting statistical 
consequence of the latter is that the polymer faces an 
entropic barrier, since the number of available configura- 
tions is least when the chain is halfway through the hole. 
In this regard it shares similarities with other entropi- 
cally controlled polymer systems, e.g. polymer trapping 
in random environments 1^-0], DNA gel electrophoresis 
1^ or reptation In these cases, the geometry of the ob- 
stacles around which the polymer must diffuse constrains 
the kinetics of the process. 

A number of recent theoretical works have shed light 
on the translocation process [^,^-14|, mostly in the pres- 
ence of a driving force. A common approach is to fo- 
cus on the dynamics of a single variable representing the 
monomer number at the pore [ p^|jT^Jl5[| . Due to its re- 
semblance to the 'reaction coordinate' for chemical pro- 
cesses, we shall refer to this parameter as the transloca- 
tion coordinate. Assuming that the segments on the two 
sides of the hole are in equilibrium leads to a force acting 



on the trapped monomer which can be derived from the 
entropic barrier mentioned before, as well as any chemi- 
cal potential differences that may provide a driving force. 
The translocation problem is thereby reduced to the es- 
cape of a 'particle' (the translocation coordinate) over a 
potential barrier. 




FIG. 1. Schematic representation of a A'^-monomer poly- 
mer in the process of translocation through a hole of size w. 
The 'translocation coordinate' s is the number of monomers 
on one side. 



Assuming Brownian dynamics for the translocation co- 
ordinate, and in the absence of a driving force, the char- 
acteristic first passage time scales as iV^, where N is 
the number of monomers. This result cannot be rec- 
onciled with the equilibration time of a polymer which 
scales as iV", with a dynamic scaling exponent a that is 
larger than 2 for Rouse dynamics of self-avoiding chains. 
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Clearly, we expect the constraint of passage through a 
hole to slow down rather than speed up the dynamics 
of the polymer, conceivably leading to an exponent even 
larger than a. Previous work recognizes this difficulty; 
for example, Ref. [T^ while using the Brownian particle 
analogy, clearly explains why it is applicable to the pre- 
asymptotic region of interest in that paper. In this work 
we consider the true asymptotic scaling of the transloca- 
tion time for large N, and find that its scaling is similar to 
the corresponding equilibration time, albeit with a larger 
prefactor. Since we cannot rely on the Brownian parti- 
cle picture in this regime, we reach this conclusion by 
numerical simulations. 

The problem is to calculate the time required for a 
polymer to move from one side of a rigid wall to the 
other through a narrow hole. This is schematically de- 
picted in Fig. with the space on both sides of the wall 
being infinite. Although frequently a driving force, such 
as an external field or chemical potential difference is 
present in the problem, we shall restrict ourselves to a 
model without external forces. For translocation to oc- 
cur, there must be two events. The first is the collective 
diffusion of the polymer to the vicinity of the pore; the 
second is its threading through the pore. For a finite 
system, or in the presence of a finite concentration, the 
first event takes a time determined by the concentration 
of polymers in solution, their diffusion constant, and the 
effective cross-section of the hole. This time is decoupled 
from the time for the second event, which is constrained 
by the passage of all monomers through the hole. Since 
we are interested only in the latter event, we shall assume 
that in the initial state the first monomer of the polymer 
chain is already threaded through the hole. To avoid the 
situation in which the polymer withdraws from the hole 
and drifts away to infinity, we add the restriction that the 
first monomer is never allowed to cross back out of the 
hole. These constraints effectively isolate the transloca- 
tion time from the time for the polymer to find the hole 

Despite the conceptual simplicity of the translocation 
problem, it has been difficult to solve analytically. Even 
the simplified case of a Gaussian polymer in a one- 
dimensional space moving past a potential barrier is non- 
trivial p6[ |. Consequently, theoretical treatments resort 
to approximations such as reducing the problem to Brow- 
nian dynamics of the translocation coordinate. As re- 
viewed in Sec. ||, the focus of this approach is the proba- 
bility density function p{s,t) that a particular monomer 
(labeled by its sequential number s along the chain) is 
located at the hole at time t. For such a constrained con- 
figuration, one can derive the entropy of the polymer if 
the segments on the two sides are in equilibrium. This 
entropy is then assumed to generate a force acting on the 
monomer, favoring its motion to one side or the other. 
Naturally, stochastic forces are also present (and are in 
fact necessary to push the chain over the entropic bar- 
rier). Assuming that the translocation coordinate obeys 
Brownian dynamics leads to a Fokker-Planck equation for 



the evolution of p(s, t). The standard Kramers' approach 
to escape over a potential barrier yields a mean translo- 
cation time that scales with the number of monomers as 
iV^, i.e. the entropic barrier does not modify the diffu- 
sive scaling. We supplement this result with a numeri- 
cal integration of the Fokker-Planck equation that yields 
the complete distribution function for transit times. The 
Rouse model for the dynamics of a phantom chain also 
predicts a time of order of N'^ for the equilibration of the 
polymer. Since such equilibration is essential to the use 
of the entropy function as the driving force, the internal 
consistency of the approach is in doubt. This is even 
more so when considering self-avoiding polymers, where 
the relaxation times are much larger than N'^. Noting 
this contradiction we proceed to numerical simulations 
in the re st o f the paper. 

In Sec. in we report on simulations of a phantom chain 
in one dimension. In this case the relaxation time of 
the chain scales as N'^, and is thus not inconsistent with 
the predictions of the Brownian particle analogy. We 
do indeed find that the probability density function for 
translocation times (once appropriately scaled) is quite 
similar to that obtained from solving the Fokker-Planck 
equation. However, the mean translocation time, while 
appearing to scale as N"^, has a much larger prefactor, 
which depends on the size of the hole. In one dimension 
self-avoiding chains arc fully stretched, and do not pro- 
vide a fair model for translocation of a coiled polymer. 

are 



Two dimensional polymers as described in Sec. IV 



ideal for studying the scaling of translocation times in 
a more complex situation. Simulations are faster and 
easier than in three dimensions, while at the same time 
the effects of self- avoidance are more pronounced. Our 
simulations for phantom chains reproduce the trends ob- 
served in one dimension, i.e. a mean translocation time 
scaling as iV^ with a larger prefactor. However, once 
self- avoidance is included the translocation times increase 
dramatically, the mean translocation time appearing to 
scale as N", with a w 2.5 which is the exponent for Rouse 
relaxation of self-avoiding chains in two dimensions. We 
find that (for the parameters of our model) translocation 
times are roughly ten times longer than typical equili- 
bration times. Thus translocation is indeed much slower 
than diffusion of the polymer, but appears to scale with 
the same exponent. 

Consequences of this observation are discussed in 
Sec. 0. The observed nontrivial scaling of translocation 
times is clear indication of the failure of the Brownian 
picture for the dynamics of the translocation coordinate. 
Instead, we suggest that the anomalous dynamics of a 
specified monomer in a chain provides a better analogy. 
Following scaling arguments used for the latter, we pre- 
dict anomalous behaviors for the translocation coordi- 
nate s{t). In particular, fluctuations in s are predicted 
to scale as f^, while in the presence of a chemical po- 
tential A/j,, the pulling velocity is predicted to behave 
as u (X A/LtiV, with ( « 0.46 and ij « —0.18 for Rouse 
dynamics in three dimensions. 
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II. BROWNIAN TRANSLOCATION 

The reduction of the translocation problem to the 
Brownian dynamics of a single coordinate was introduced 
in Ref . [0 , and further explored in Ref . . Here we re- 
view the main features of this approximation and its con- 
sequences. Consider a polymer moving through a pore 
in a membrane, where the hole is so narrow that only a 
single strand of polymer can pass through. (Thus, the 
parameter describing the width of the hole in Fig. ^ is 
w = 1.) The progress of the polymer can be tracked by 
following the number s of the monomer which is located 
in the hole at a particular time, as depicted in Fig. 0. Let 
us denote the probability of monomer s being in the hole 
at time t by p{s,t). As the monomer s moves forward 
or backward through the hole a distance a (of order of 
the typical separation between monomers) the relevant 
monomer number increases or decreases by unity. Treat- 
ing s as a continuous variable, we can write a continuity 
equation for the probability as 



dp dj_ 
dt ds 



0, 



(1) 



where j(s, t) is the probability current. 

The central difhculty is to find an appropriate expres- 
sion for j which correctly reflects the correlated motion of 
the whole polymer. If the progress of the polymer is suf- 
ficiently slow for the segments on the two sides to come 
to equilibrium, the monomer at the hole experiences a 
mean force that can be obtained from the variations of 
a constrained free energy F(s). How such a force can be 
used to deduce the dynamics of the monomer label is not 
clear. The analogy to Brownian motion suggests that the 
rate of change of s is related to the force by a mobility 
fj,. Since the polymer fluctuates back and forth between 
the two sides, there must also be a stochastic element 
which can be represented by a random force. If there are 
no correlations in this force at different times (as in the 
standard Langcvin formulation for a Brownian particle) , 
there is a current that depends on the local probability 
density as 



J 



dp p OF 



(2) 



where the diffusion parameter D is proportional to the 
variance of the stochastic force |18|. As in standard 
Brownian processes, the above equation assumes that the 
mobility is related to the temperature T hy fi — D/ksT, 
where ks is the Boltzmann constant. No similar restric- 
tion is made in Ref. [ p^ , which obtains the diffusion term 
from considerations of symmetry and locality. However, 
as elaborated in Sec. ^ the assumption of locality need 
not be valid in this case, since the true dynamics of s must 
reflect the collective behavior of the whole polymer. 

Calculating the restricted free energy F{s) is reason- 
ably straightforward, and equivalent to finding the num- 
ber of possible configurations of a polymer attached at 



one point to an impenetrable barrier. The exact solution 
is known for the case of one-dimensional discrete random 
walks with fixed step-length (see, e.g., Ref. (l^): It can 
be shown that (in the large N limit) the number of iV- 
step walks which start at a boundary and never return 
to it is yj2/TTN ■ 2^ . Thus, the number of configurations 
with s monomers on the right and N — s monomers on 
the left has the s-dependence A/[{N — s)s]'' , where A is 



independent of s and 7 = 
part of the free energy as 



1/2, giving the s-dependent 



F = 7/cBTln[(7V-s)s]. 



(3) 



To this result we add the conditions that the first 
monomer can never be withdrawn from the hole, and 
that after the JVth monomer crosses the wall the poly- 
mer will no longer return to it. Figure ^ depicts the 
resulting free energy for the case of iV = 1000. The two 
conditions are shown by two vertical lines on the sides 
of the graph: the line on the left (infinite barrier) sig- 
nifies our assumption that the first monomer can never 
cross back through the hole, while the vertical line on 
the right (—00) represents the escape of a polymer which 
has crossed the barrier. The numbers of configurations of 
polymers in higher space dimensions d, whether phantom 
or self-avoiding (SA), cannot be calculated exactly. How- 
ever, it is known that they have the same dependence on 
the polymer length with 7 in Eq. (^) replaced by 
a number which depends on rf (7 = d/2 for phantom 
chains). Thus the s dependent part of the free energy 
will have a different prefactor, but the logarithmic de- 
pendence will remain unchanged. 




FIG. 2. The entropic potential barrier found in the single 
variable equivalent to the polymer translocation problem. 

By inserting Eq. (^) into Eqs. (|l]) and (||), we reduce 
polymer translocation to a standard single particle prob- 
lem of escape over a potential barrier. It is interesting to 
note that if we rescale the variables according to s — > sN, 
t tD/N'^, all iV- and D-dependence is eliminated from 
the equation, resulting in 
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dp _&^ d_ / l-2s \ 

The solution to this dimensionless equation can be con- 
verted back to real time by multiplying the t axis by 
N"^ /D. Thus, under the assumptions listed above, the 
escape time of a particle, and thus the translocation time 
of a polymer are proportional to N"^ /D. Note, that this 
conclusion is independent of the value of the parameter 
7, and remains valid even for a SA chain in which 7 has 
a different value. 

The problem of escape from a deep well in one dimen- 
sion was considered by Kramers Assuming that 
the escape rate is slow, i.e. if at every moment in time 
the probability distribution of a particle in the well can 
be represented by an equilibrium (Boltzmann) weight, 
Kramers' method enables an analytic calculation of the 
mean escape time. Applying Kramers' formula to the log- 
arithmic potential of the problem one finds (for 
7 — 1/2) that the mean escape time t is {tt'^ /16)N'^ / D. 
The distribution of escape times in Kramers' formula is 
by construction a simple exponential. Equation ^ can 
also be solved numerically, by placing a delta function at 
the left edge of the potential depicted in Fig. || at i = 
and integrating in time. There are some differences be- 
tween the numerical solution to Eq. ^, and the Kramers' 
solution, e.g. the former decays to zero for t ^ due to 
the time it takes for the delta function to diffuse out of 
the well. However, the mean escape time for the numeri- 
cal solution is r w 0.6N^/D which almost coincides with 
the approximate Kramers' result. The distribution of es- 
cape times calculated from integrating the Fokkcr-Planck 
equation is shown by the dashed line in Fig. ^. 




FIG. 3. Probability distribution of escape times measured 
in reduced, dimensionless, units. The dashed curve cor- 
responds to the logarithmic potential depicted in Fig. ^ 
while the solid line is in the absence of the entropic bar- 
rier. The curves were obtained by numerical integration of 
the Fokker-Planck equation. 

If the potential barrier is absent altogether, then the 
distribution of escape times can again be calculated (e.g. 



by numerical integration), and is depicted by the solid 
line in Fig. ^. The mean escape time in this case is 
r « 0.5 (in reduced units of N'^/D). There is only a 20% 
difference between the mean escape times of the prob- 
lems with and without the logarithmic potential barrier. 
More strikingly, there is little difference in the distribu- 
tion of times with or without a barrier. (This barrierless 
version of the problem can be thought of as describing 
the adsorption of a particle starting at a unit distance 
from a sink. From a polymer perspective, it corresponds 
to passage of a polymer through a ring - although the 
polymer must pass through a constricted space, its free 
ends can have any possible configuration.) In view of the 
minute effect of the entropic potential, the notion of 'es- 
cape over a barrier' does not provide a particularly useful 
analogy. 

A central assumption in the reduction of the polymer 
problem to a single coordinate is that translocation is 
slow. Specifically, it should be slow enough that the poly- 
mer segments on the two sides of the membrane are in 
equilibrium at every value of s. We can check for the self- 
consistency of this assumption: The equilibration time of 
a free polymer can be estimated |^ as the time required 
for it to diffuse a distance equal to its own radius of 
gyration Rg. Under Rouse dynamics (which ignores hy- 
drodynamic effects), the diffusivity of the center of mass 
of an iV-monomer polymer is reduced to D/N, where D 
is the diffusion constant of a single monomer, resulting in 
an equilibration time of order of R^N/D. For phantom 
polymers R^ ~ N, and, consequently the equilibration 
time is of the same order (N'^/D) as the mean passage 
time obtained with Brownian translocation dynamics. 
Thus the equilibrium assumption is marginally (in)valid 
in this case. On the other hand, for self-avoiding poly- 
mers Rg ~ iV", where v = 0.75 and 0.59 for dimen- 
sions d = 2 and 3, respectively. The resulting relaxation 
times (^ N^^'^'^/D) are now longer than the transloca- 
tion times predicted (~ N'^ / D) by Brownian dynamics, 
and consequently the approximations involved are not 
self-consistent. In the following sections we try to gain 
further insights into the problem by numerical simula- 
tions of the translocation of a polymer in one and two 
dimensions. 



III. SIMULATIONS IN ONE-DIMENSION 

We begin by examining the translocation of a one di- 
mensional phantom polymer via a Monte Carlo simula- 
tion with Rouse- like dynamics as follows. Our model 
consists of a chain of N 'atoms' placed on the sites of 
a one-dimensional lattice. No excluded volume interac- 
tions are present, and the spatial distance between two 
neighboring atoms (along the sequence of the chain) can 
be 0, 1 or 2 lattice spacings, i.e. the "bond" between ad- 
jacent atoms has a maximal length of 2. This represents 
a trivial implementation of the fluctuating-bond method 
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| p2[ . Initially, the first atom of the chain is placed on, say, 
the right of the membrane, while all other atoms are on 
the left. (We assume that the membrane is located be- 
tween the coordinate x = and a; = 1 and, consequently, 
the first atom is initially placed at x = 1 while the rest of 
the atoms are at a; < 0.) During the simulation the first 
atom of the chain is never allowed to move to the left of 
the membrane. The "width" w of the hole is adjusted 
by changing the maximal number of bonds allowed to be 
simultaneously present at the hole. An elementary move 
consists of randomly picking an atom and attempting to 
move it one lattice step in a randomly selected direction. 
If the new configuration does not violate any of the re- 
strictions of the model it is accepted. N elementary atom 
move attempts are defined as one Monte Carlo time unit. 




12 3 

l./T 



FIG. 4. Probability density of translocation times for 
one-dimensional phantom polymers. Solid, dotted and 
dashed lines correspond to polymer lengths — 32, 64 and 
128, respectively, and were extracted from 10,000 simulations, 
each. The times have been normalized by their respective 
mean translocation time t{N) in each case. The continuous 
curve corresponds to the solution of the Fokker-Planck equa- 
tion in the single-particle approximation. 

Each simulation is terminated when all monomers are 
on one side. Figure ^ depicts the distribution of translo- 
cation times of such phantom chains measured for sev- 
eral chain lengths, and for unit width of the hole. If we 
normalize the time for each N by the mean transloca- 
tion time T for that length, we observe that the resulting 
curves are quite similar and closely resemble the theoret- 
ical curve obtained assuming a Brownian translocation 
coordinate. While the similarity may appear to support 
this picture, we should note that since the distribution is 
constrained to vanish at both short and long times, qual- 
itative similarities are dubious. Moreover, the absolute 
values of the mean translocation time r, as depicted in 
Fig. 1^, are significantly larger than the estimates from 
Brownian dynamics. This log-log plot indicates that the 
apparent exponent is somewhat larger than 2 for small TV, 
and only gradually approaches the scaling form r ~ iV^. 
However, the prefactor of the power law for N = 256 



is roughly two orders of magnitude larger than expected 
for a Brownian translocation coordinate. Such discrep- 
ancy should not be surprising - the translocation time 
predicted by this model are similar to the time required 
for a polymer to diffuse its own radius of gyration. It is 
reasonable to expect that passage through a narrow hole 
should be slower than diffusion without a wall. 
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FIG. 5. Logarithmic plot of the mean translocation time r 
as a function of chain length A'^ for a one dimensional phantom 
polymer. Each data point represents an average over 10,000 
realizations. The solid line has slope 2. 



Not surprisingly, translocation times are strongly af- 
fected by the width w of the hole, as indicated in Fig. |^. 
When w reaches the size Rg of the chain (~ Vn) the 
translocation time should become independent ofw. This 
is supported by the saturation of the rescaled times t/N"^ 
at roughly the same value of w/Rg ~ w/^/N in this fig- 
ure. For small holes, the translocation times are strongly 
dependent on the hole size, and show some indications of 
collapse onto a universal curve, but the N'^ dependence 
is not as clear as for the wider holes. 
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FIG. 6. Dependence of the mean translocation time r on 
the width of the hole w, for = 3, 4, 6, 8, . . . 91, 128, and 
181. The times have been normalized by A'^'^ to focus on the 
behavior of the prefactor. Each curve corresponds to fixed A^, 
and is obtained by averaging over 1,000 cases. 



lie in the aperture at a time, this condition is equivalent 
to the complete translocation of the polymer through the 
hole. As in the one-dimensional case, the first monomer 
is not allowed to move to the left of the hole. 



IV. SIMULATIONS IN TWO DIMENSIONS 

Excluded volume effects drastically modify the shape 
and properties of the one-dimensional phantom polymer 
considered in the previous section. The chain becomes 
stretched, and its dynamics are then limited by reptation 
1^ . To study the effects of self-avoidance on translocation 
in the coiled state, higher dimensional simulations are 
necessary. Two-dimensional polymers are ideally suited 
to this purpose for the dual reasons that excluded volume 
effects are more apparent, while computation times are 
shorter than in the three-dimensional case. 

As in the 1-d system, we employ a fluctuating-bond 
model for simulations p^ , implementing Rouse- like dy- 
namics for linear polymers of several different lengths N. 
Simulations were performed both with and without ex- 
cluded volume constraints. In this model, the monomers 
of the polymer lie on a 2-d square lattice. Random mo- 
tion is simulated through a series of elementary moves of 
single monomers. In each move, a monomer is selected 
randomly and then moved a single lattice unit in one of 
the +x, —X, +y, or —y directions. If the move violates 
any of several constraints, it is rejected. For phantom 
chains, bonds have a maximal allowed length of vTO lat- 
tice units. For excluded volume chains, in addition, the 
distance between any two monomers is constrained to 
be at least 2 lattice units. For excluded volume chains, 
these constraints also prevent the chain from crossing it- 
self . The wall has a thickness of 3 units and the hole 
has a width of 2 lattice units. At these sizes, only one 
monomer may be in the hole at a time, but the hole is 
large enough that translocation can occur with the given 
move set. Each N elementary atom move attempts are 
defined as one Monte Carlo time unit. 

The simulation begins by placing the first monomer 
at the hole, while the remaining A'^ — 1 monomers are 
in a random conformation on the left side of the wall. 
(To generate the initial random configuration, the chain 
is first allowed to fluctuate subject to the constraints of 
impermeability of the wall and fixed location of the first 
monomer. The fluctuation time is 2QN'^, which should be 
sufficient to randomize the initial condition for the pur- 
poses of our simulations, at the length scales employed. 
In any case, the time it takes for translocation is many 
times longer than the Rouse relaxation time, and any 
initial condition effects will be transients.) Once the ini- 
tial configuration has been established, the polymer is 
allowed to move in accordance with the restrictions of 
the model. We measure the time between the begin- 
ning of the translocation, and the moment when the last 
monomer enters the hole. Because only one monomer can 



le+07 
le+06 
le+05 
10000 
^ 1000 
100 
10 
1 



O Excluded Volume Chain 
O Phantom Chain 



16 32 64 128 



FIG. 7. Logarithmic plot of the mean translocation time r 
as a function of chain length A*', for two-dimensional phan- 
tom and excluded volume polymers. Error bars indicate the 
standard deviation over runs. 



Translocation times were calculated for a number of 
different chain lengths A^, with several thousand runs at 
each length. (The number of runs decreases with increas- 
ing A", due to CPU limitations.) Theses results are shown 
for both phantom and self-avoiding chains in Fig.|^. Er- 
ror bars indicate the standard deviation of translocation 
times over runs. 
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FIG. 8. The effective power law exponents for translocation 
and diff'usion times as a function of length N for 2-d polymers. 
For both translocation and free diffusion, self-avoiding chains 
exhibit a dynamic exponent that approaches 2.5, equal to 
1 + 2v given the d = 2 swelling exponent of i/ = 3/4. For 
phantom chains, the exponent is noticeably smaller, and ap- 
proaches 2. In both the self-avoiding and phantom cases, the 
asymptotic exponent for translocation time is nearly equal to 
that for Rouse relaxation. 
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We then attempted to fit the results to a dynamic 
scaling form, with the effective exponent depicted in 
Fig. ||. The points are a plot of the local exponent a, 
assuming a power law r ~ iV", as determined from two 
neighboring polymer lengths A'^i and N2 via the formula 
a(V^VT7V2) = log[r(iV2)/T(Ari)]/log(iV2/A^i). The expo- 
nent a for the excluded volume translocation approaches 
an asymptotic limit which is definitely larger than 2. This 
clearly indicates that excluded volume effects are incon- 
sistent with Brownian dynamics for the translocation co- 
ordinate. The data for the translocation times of a phan- 
tom polymer (in which the self-avoiding constraint of the 
bond-fluctuation model is removed) are also included in 
this figure. The power law exponent in this case asymp- 
totically approaches a value close to 2, in agreement with 
the one dimensional results, and consistent with Brown- 
ian dynamics. At each iV, the uncertainty in the average 
translocation time (equal to the standard deviation of 
t{N) over runs divided by the square root of the number 
of runs) has been translated into an uncertainty in the 
local power law. This uncertainty, as depicted by the er- 
ror bars, is quite small, indicating that our conclusion is 
not due to statistical fluctuations. 

The observed scaling laws for self-avoiding and phan- 
tom polymers in fact agree with the exponents expected 
theoretically for Rouse relaxation in the absence of a wall. 
As stated previously, translocation across the barrier re- 
quires that, minimally, the chain diffuses a distance equal 
to its radius of gyration. The time for such diffusion is 
Tfl - RlN/D - For a self-avoiding walk in 2 

dimensions, v — 3/4 and ~ N"^-^ , while for a phan- 
tom chain, v is replaced by 1/2, leading to the relation- 
ship Tn ^ N"^. For comparison. Fig. ^ also shows the 
derived power law exponents from simulations of simple 
diffusion in the absence of a wall. For an excluded vol- 
ume chain, we approach the scaling r ~ A^^'^, while for 
the phantom chain, r ^ A^^. The latter results were 
generated as follows: We started with the same type 
of random initial configurations (generated by annealing 
near an impermeable wall) as for the translocation case, 
and again imposed the constraint that the first monomer 
cannot pass to the left side of an imaginary wall. The 
other monomers, however, are allowed to diffusive with- 
out feeling this wall. The simulation was stopped when 
all monomers moved to the right of the imaginary wall. 

At least in these examples, we find that the scaling of 
translocation times is the same as that for equilibration 
of the polymer in the absence of a wall. However, the 
constraint of passing through the hole must clearly slow 
down the dynamics of the polymer compared to the case 
of free diffusion over a similar distance. This slowdown 
must then be reflected in an overall prefactor which deter- 
mines how much slower translocation is relative to pure 
diffusion. The data that address this issue are displayed 
in Fig. ^ The ratio of crossing times for translocation 
compared to pure diffusion (as described in the previous 
paragraph) are plotted for both self-avoiding and phan- 



tom chains. For self-avoiding polymers, the hole slows 
down the chain by a factor of about 5; while for phantom 
chains the ratio is about 13. These numbers are roughly 
A^-independent for A^ > 10 (for excluded volume chains) 
and A^ > 50 (for phantom chains). It is interesting to 
note that translocation slows down the phantom chain 
more than it does the excluded volume chain. However, 
these ratios should not be taken too seriously, since for 
more realistic systems such as the translocation of DNA 
molecules through a cell membrane, the details of the 
shape and interaction forces at the pore play a signifi- 
cant role lH,pT 



• Time Ratio for Phantom Chains 

O Time Ratio for Excluded Volume Chains 



FIG. 9. Comparison of simulation results for diffusion 
through a hole that can fit one monomer versus free diffu- 
sion. Translocation times through the hole are larger than 
the free diffusion for both phantom and self-avoiding chains. 
In both cases, the ratio of these times levels off to a constant. 



V. DISCUSSION 

The central result of this paper is the nontrivial scaling 
of the mean translocation time r for unforced passage of 
a polymer through a hole with its length A^ . For the diffu- 
sive dynamics of a self- avoiding chain in d = 2, numerical 
simulations indicate r ~ A^^'^. By extension, we expect 
T ~ A^" with a = l + 2i> for diffusive dynamics. Possibly 
even more generally for other types of dynamics a = 1^2, 
where the swelling exponent v relates the radius of gy- 



ration of the polymer to its length by Rg 



N", while 



its relaxation time scales as Tr ~ Rg with the dynamic 
exponent z. 

The natural parameter for following the progress of the 
translocation process is the label s(t) for the monomer in 
the hole at time t (see Fig. ||). A commonly used anal- 
ogy is that this 'reaction coordinate' undergoes stochas- 
tic motion, as in a Brownian particle in a force field 
[]lO| , ^ 13 1 . For unforced translocation (with or without 
inclusion of an entropic barrier), such approaches lead 
to T ~ A^^, which is inconsistent with our numerical re- 
sults, except for the case of diffusing phantom polymers. 
Such scaling is also at odds with the expectation that 
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constraining the polymer to pass through a hole should 
slow down its dynamics compared to free motion over 
a similar distance. The numerically observed scaling is 
thus indicative of anomalous dynamics of the transloca- 
tion coordinate s(t) — we propose instead an alternative 
analogy that incorporates this observation. 

The dynamics of the polymer is the result of coopera- 
tive motions of its many monomers. When described by 
just a single variable, say the translocation coordinate, 
the effect of the other degrees of freedom is to exert cor- 
related forces leading to anomalous dynamics. This is 
best illustrated by the well studied case of the dynamics 
of a single monomer in a polymer: Consider the posi- 
tion f{t) of a particular atom in the background of all 
the other monomers. The dynamics of r{t) has indeed 
been studied by numerical, analytical, and experimental 
methods |^ ; its anomalous features are easily captured 
by the following scaling argument: For short times, we 
expect that the squared change in the position has a scal- 
ing form (Ar^(i)) ^ i^'', with no dependence on N since 
the monomer does not yet feel any effects from the finite 
size of the polymer. At a time of the order of the relax- 
ation time Tr for the whole chain, the monomer should 
have moved by a distance of the order of its radius of gy- 
ration, giving i?g ^ T^'>. Since ~ Rg, we immediately 
obtain the exponent C. — l/z describing the anomalous 
fluctuations of the specified monomer at short times. 

We now adapt a similar scaling argument to describe 
the squared change in the translocation coordinate, as- 
suming (As'^it)) - at short times. This behav- 
ior should saturate when s becomes of the order of the 
chain length in a time t ^ Rg ^ N^^ (assuming that 
translocation times always scale in the same way as equi- 
libration times). Substituting this in the former equa- 
tion allows us to identify the exponent C = l/{zi') for 
anomalous dynamics of the translocation coordinate. For 
the case of diffusive polymer dynamics, we thus obtain 
C = 1/(1 -I- 2v) resulting in C = 2/5 and C « 0.46 for self- 
avoiding chains in two and three dimensions respectively; 
i.e. in this case the fluctuations are subdiffusive. If we 
naively use the Zimm exponent z — d describing relax- 
ation of polymers in hydrodynamic flows in d dimensions, 
we obtain l{d = 2) = 2/3, while C{d = 3) w 0.56, i.e. 
the fluctuations are predicted to be superdijfusive in this 
case. The usual origin of the speed up of polymer dy- 
namics in a fluid is attributable to the velocity flow field 
set up by the other monomers. It is doubtful that using 
the bulk Zimm exponent remains valid for flows which 
must vanish in the vicinity of the wall. Thus the above 
prediction of superdiffusive behavior should not be taken 
seriously prior to a proper analysis of the hydrodynamic 
correlations in the vicinity of the wall. 

It is experimentally hard to directly probe the motion 
of the translocation coordinate. The quantity that is eas- 
ily measured in experiments Q is the distribution of the 
translocation times in the presence of a force introduced 
via a chemical potential difference A/j, for monomers on 



the two sides of the wall. While our results so far were 
in the absence of such a driving force, anomalous dy- 
namics has consequences for the length dependence of 
the forced velocity. Let us first recall the arguments for 
the drift velocity f of a polymer in a force F: Scaling 
considerations suggest v{F) ^ {Rg/Tr)(f>{FRg/kBT) ^ 
7V''(i-^)(/.(FiV''), where is a scaling function depend- 
ing on the ratio of two quantities having dimensions of 
energy. The proportionality of the velocity to the force 
requires a linear scaling function, leading to a mobility 
w/F - N-i^i^-^), For Rouse and Zimm dynamics this 
leads to the well-known scalings of mobility as 1/A^ and 
1/Rg~^, respectively. Similarly, for the pulling velocity 
M = s of the translocation coordinate, scaling suggests 
u ~ N/T(l)iAfiN/kBT) - N^-^^Afi. Only for the case 
of diffusive dynamics of a phantom chain this velocity 
is independent of N. The anomalous slowdown due to 
Rouse dynamics leads to a mobility that scales as N~^^'^ 
in d — 2, and N^'^-^^ in = 3. By contrast, hydrody- 
namic speeding up with Zimm dynamic exponents leads 
to a mobility that grows as N^/"^ and iV^-^^ in two and 
three dimensions respectively. Once more, the latter re- 
sults are not to be taken seriously without full hydrody- 
namic calculations in the presence of a wall. Experiments 
so far Q do not indicate anomalous scaling, but the size 
range may not be sufficient to detect the rather small 
exponent. 

We conclude by listing two other avenues of poten- 
tial exploration. The first is to note that the distribu- 
tion of translocation times should also be modified by 
the anomalous dynamics, potentially to include power 
law tails that are distinct from the exponential tails in 
Fig. H p^ , p5t . Secondly, it may be possible to construct 
an experimental system close to our two dimensional sim- 
ulations, using vibrated granular chains in a variation on 
the set up used in Ref. |26|| , with a chain that is threaded 
through a hole in a wall. 
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